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ABSTRACT 

Side discharge of power plant thermal effluents into natural rivers 
cause the formation of thermal plumes, the shape and orientation 
of which depend upon both the discharge conditions and the ambient 
flow characteristics. When the turbulence level of the ambient flow 
is high, rapid vertical mixing is achieved between the effluent 
discharge and the river flow, and the thermal plume becomes shore- 
attached with the maximum temperatures within the plume occurring 
along the near shore. In order to determine the isotherm patterns 
within shore-attached thermal plumes, a solution of the depth-inte- 
grated steady state convection-diffusion equation is obtained. The 
solution is cast in terms of space-cumulative discharge coordinates, 
thus allowing the inclusion of the effects of cross-channel flow 
induced by channel curvature in the model. Comparison of the model 
results with the thermal plume data from power plant sites along 
the Missouri River indicates good agreement. 


INTRODUCTION 

In order to meet the increasing electrical energy requirements of 
the nation, steam-electric power plants of large capacities are 
being installed or being planned for installation by a number of 
utilities. However, the low thermal efficiencies of these plants 
necessitate the rejection of large amounts of waste heat at the 
generating sites. In situations where it can be successfully 
demonstrated that discharge of the waste heat into a natural body 
of water will not cause any harmful effects on the ecology of the 
receiving aquatic environment, open-cycle cooling using river, 
lake, or sea water is the most desirable waste-heat-removal method. 
Closed-cycle cooling systems using cooling towers, cooling ponds, 
or spray canals involve enormous capital investments for construc- 
tion; these can be up to 15 times that for an open-cycle system. 
They may also require a considerable portion of the plant output 
power for operating the cooling systems, thereby imposing penalties 
on plant efficiency. Closed-cycle systems also require extensive 
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continuing maintenance. Finally, there is consumptive water loss 
due to evaporation from these systems. Open-cycle systems, on the 
other hand, are inexpensive to construct, place minimal demands on 
the plant output, and evaporative water losses are much less. Thus, 
lower investment and maintenance costs, improved plant efficiencies, 
and less water loss all combine to make once-through cooling the 
most attractive method by far, from an economic point of view, for 
dissipating power plant waste heat. 

The effects of the thermal discharge on a natural stream, and the 
physical processes involved in heat transport within the flow system 
are normally analyzed by separating the affected area into two re- 
gions; the far-field and the near-field. The determination of far- 
field temperature effects in a river is important in assessing the 
cumulative effects of several power plants distributed along the 
course of the river. On the other hand, the determination of near- 
field conditions is important in evaluating the possible effects of 
the thermal discharges on the aquatic biota and drifting organisms. 
In the vicinity of an outfall, the heated water forms a so-called 
thermal plume until the combined actions of jet mixing, buoyant and 
convective spreading, and turbulent mixing eventually lead to com- 
plete dilution of the heated water with the ambient flow. The ther- 
mal standards of the various regulatory agencies regarding thermal 
discharges in natural waterways usually specify certain limitations 
on the size of and temperature rise in thermal plumes. It is there- 
fore important to have reliable methods for predicting configura- 
tions of and temperature patterns within thermal plumes. 

Several computational methods are available for predicting thermal 
plume patterns resulting from power plant discharges. Recently 
Jirka, Abraham and Harleman (4) have published a detailed review 
of these. Most of these techniques either assume or predict bell- 
shaped profiles for velocity and temperature along horizontal nor- 
mals to the plume axis. However, an examination of measured ther- 
mal plume patterns near power stations at Omaha, Fort Calhoun, and 
Brownville (Cooper) on the Missouri River indicates that the tem- 
perature profiles are not beel-shaped, but are attached to the 
shores, with the maximum temperatures in the plumes occurring near 
or along the shore. None of the currently available thermal plume 
models predict plume patterns with shore attachment; this deficiency 
is highlighted by Jirka, Abraham and Harleman (4). 

The available the*rmal plume models fail to properly represent the 
shore-attached plume patterns because the effects of ambient flow 
characteristics such as bottom effects, variation of channel shape, 
non-uniform and varying flow conditions across the channel, and 
ambient turbulence on the mixing phenomenon are not considered in 
their formulation. If the turbulence level of the ambient river 
flow is high, (e.g., the Missouri River), intense mixing between 
the thermal jet and the ambient flow will be induced, leading to 
rapid vertical mixing between the two flows. Thereafter, the jet 
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is transported downstream as a coflowing stream, and the mixing is 
mainly due to turbulent diffusion and transverse convective spread- 
ing associated with bend-generated secondary circulation. Weil and 
Fischer (13) have attempted to relate the effects of stream tur- 
bulence to plume behaviour. However, their investigation was re- 
stricted to zero-momentum jets discharging at the centerline of the 
ambient flow. Yeh (14) conducted a series of laboratory experiments 
to investigate the transverse mixing process for a buoyant effluent 
in both near and far field regions of an open channel flow. The 
results showed that the mixing process in the near field region was 
closely related to the development and decay of buoyancy-generated 
secondary flow. However in the far field region the transverse mix- 
ing was found to be mainly due to turbulent diffusion. 

The present study is concerned with the development of a predictive 
model, based on diffusion concepts, for determining transverse tem- 
perature profiles in shore-attached, vertically-mixed thermal plumes 
in rivers resulting from side discharges of thermal effluents. The 
details of the analytical development of the model, and its applica- 
tion to field cases are presented in the following sections. 


FORMULATION OF DIFFUSION PROBLEM 

The formulation of the predictive model for determining the temper- 
ature distributions in shore-attached plumes is based on a model of 
mixing of two rivers outlined by Sayre (9) and the transformation 
of the two-dimensional depth-integrated steady-state convection- 
diffusion equation into a simple diffusion equation as presented 
by Yotsukura and Sayre (17). The solution of the transformed equa- 
tion is cast in terms of space-cumulative discharge coordinates, 

(x, p(x,z)), when x and z are the streamwise and cross-channel coor- 
dinates, and p(x,z) is the normalized cumulative discharge defined 
by 



( 1 ) 


where Q„ = total river discharge; q = cumulative discharge passing 
through^the river between z = 0 and^z at section x; and q = total 
discharge per unit width. The use of the space-cumulative discharge 
coordinate system which automatically follows the transverse shifts 
and meander ings of the flow, was first introduced by Yotsukura and 
Cobb (15). A procedure for synthesizing the normalized cumulative 
discharge at cross sections of a river for which the shape is known 
is outlined in the next section. 
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The origin of coordinates, x=0,p«0, is located on the discharge 
bank of the river, just far enough downstream from the outfall so 
that at the origin initial mixing has taken place, and the tempera- 
ture rise of the jet is 0-. The fraction of the total river dis- 
charge occupied by the jet at the origin is given by 


P 





(2) ^ 


where a = ®o/®t “ initial dilution factor? 0^ = initial excess 
temper ature°of‘'the thermal discharge? and = volumetric rate of 
effluent discharge. 

« 

Let the variable ? represent the displacement from the origin in the 
p direction at x = 0, as shown in Fig. 1. Then the input transverse 
temperature distribution at x = 0 is given by 





0 < ? ^ p'\ 


0, P < C < 1 


(3) 


Neglecting for the moment reflections from the banks, the normalized 
transverse temperature distribution function for the case of a con- 
tinuous vertical line source, extending over the depth of flow and 
concentrated at (x = 0, p = c) is f„{p- 5 ? x) which is a Gaussian pro- 
bability distribution function, given by 


^ ^‘1 Pi = P~^ 


where s = p^/a is the standardized normal variable? and the standard 
deviation in tSe p domain is 


Op = V2D X ' 


(5) 
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In Eq.(5), D = overall transverse diffusion factor; and x' = x - x 
where x » longitudinal coordinate of the virtual source for any ° 
subreacn determined from 









(6) 


] 


where ® distance from the origin to the upstream end of j-th 

subreacnT'^and D. = value of D for the j-th subreach. For the first 
subteach, x^j =^0. 

The desired output temperature distribution for an input distribu- 
tion 0) is given by the convolution integral 


-P=l 

AT(p,x) - fj^(p-?;x) 0j(^,O) 


(7) 


Using Eq. (3), Eq. (7) can be written 

AT(p,x) = fj^(p-5;x) d? 


( 8 ) 


Taking into account now reflections from the banks, the final solu 
tion can be obtained from Eq. (8) as 


AT(p,x) = 


R «p 
(- 

n=l ' R "P 


H- I { F (2n_±_i£±Zl) - p (2£L^i£±£l) 

, L D On V o- 


(9) 
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where Fj^( ) is the standardized cumulative normal distribution 
function corresponding to the probability density function f-( ). 

The terms in the infinite series representing the bank reflections 
are negligible if Dx< 0.08. 

The value of the initial dilution factor, a, introduced in Eq. (2) 
is mainly a function of the design and orientation of the discharge 
outfall structure, the initial jet or plume characteristics, and 
the ambient flow properties. If the thermal effluent is discharged 
almost parallel to the river flow, and ©^is the temperature excess 
measured at the discharge canal outlet, the initial dilution factor, 
a, can have a value of one. For other situations, values of the 
initial dilution factor between about 2 and 4 can be used (3,7), 
until better criteria are established. 

Maximum Centerline Temperature 


The maximum centerline temperature, 
occurs along the near shore (p = 0). 
Eq. (9) is 


AT , for a shore-attached plume 

The solution for AT^, from 

c 


AT = AT(0,X) 
c 



+ 2 


I 

n=l 





2n+P 
o 



. 2n-P 


)} 



(10) 


The variations of values of P equal to 0.05, 

0.10, and 0.20 are shown in°Fig. 2. The limiting case, P = 0, 
corresponds to a source concentrated at the origin. The solution 
corresponding to Eq. (9) for this limiting case is 


AT (p,x) 


200— [fRtP) + ^ { fR^2n-p) + fR(2n+p)}) 

Qr n=l 


( 11 ) 


Substituting for the probability density function fp(p) by 


“'7^ exp(-^ , 


Op>/2x 


(12) 
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the solution for maximum temperatures along the shore; (p » 0), 
corresponding to Eq. (10) for this case is 


ATq = AT(0,x) 




(1 + 2 


7 exp (- 25 ?^)). 

n=l p 


(13) 


Eq. (13) is represented by the curve labelled concentrated source 
in Fig. 2. 

Diffusion Factor, D 

The diffusion factor, D, in Eq. (5) can be_related to the dimension 
less tranverse mixing coefficient, a= E /du*, by the equation 

Z 


D 


a 


d u» 


1 

( ^ d^ u dp) 
0 


(14) 


if it is assumed that the overall transverse mixing coefficient E 
does not vary across the channel. In Eq. (14) d = A/B = average, 
depth of flow section; A = area of river flow section; u*= (gdS)"^' = 
shear velocity; g =■ acceleration due to gravity; S = slope of energy 
gradient; d = local flow depth; and u =* q/d =* local depth-averaged 
velocity of flow. Using the Manning formula for the river flow rate. 



1-49 

n 


B d (d) ^ 



(15) 


and approximating 


\ d^ u dp ''' (3^) ^ u = d q, 
•*0 


(16) 


where q = d u = width-averaged value of q, the relation for the 
diffusion factor can also be written 


D 


a nVg 
1*49 



(17) 
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where n = Manning's coefficient. Values of D determined by Eg. (17) 
are less accurate compared to those by Eg. (14), mainly because of 
the approximation used in Eg. (16). 

The transverse mixing coefficient, E or a, for a natural river 
depends«>upon the characteristics of Ihe river channel and the 
ambient flow. However, there are as yet no very reliable predic- 
tive relationships for evaluating its value. The most reliable 
method for determining the values of E for natural channels is by 
direct experimentation in the field usJng tracers. Experiments in 
undistorted hydraulic models can also be used to determine E ; 
vertically distorted hydraulic models do not correctly simulate 
the transverse mixing processes. 

In laboratory flumes several investigators have found val^ues of 
a ranging from about 0.1 for B/d r 5 to about 0,2 for B/d - 60. 

In reasonably straight uniform reaches of natural channels some- 
what larger values of a have been found. This increase is pro- 
bably due to channel irregularities and weak secondary flows rather 
than to any scale effect. In sinuous channels, due to bend-induced 
secondary flow, avalues are typically much larger. As shown in 
Fig. 3, average avalues ranging from 0.5 to 2.5 for bends in 
flumes, and from 0.6 to 10 for bends in the Missouri River, have 
been found. The functional form of the relationships shown on 
Fig. 3 was derived by Fischer (2) in a transverse convective dis- 
persion analysis which utilized Rozovskii's (9) radial velocity 
distribution function for fully-developed secondary flow in the 
central portion of an idealized curved channel. The gap between 
the laboratory flume data and the Missouri River data in Fig. 3 
indicates that there are still certain factors that are not 
accounted for in evaluating the transverse mixing coefficient. 

The eguation for the curve passing through the Missouri River 
data is 



(18) 


where u = Qg/A = average river flow velocity; and R = radius of 
curvature or river bend. The value of the numerical coefficient 
in Eg. (18) may be different for other rivers. 


SYNTHESIS OF TRANSVERSE FLOW DISTRIBUTION 

The cumulative discharge, q , and the normalized cumulative 
discharge, p = ^ fiver section at a given flow discharge 

can be synthesized rrom the transverse depth profile of the cross- 
section, as outlined by Sayre (10). The method is based on the 
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hypothesis that the velocity distribution in natural streams is a 
function of the transverse depth profile, and uses the relation 


_SU 


q 


^o 



(19) 


where d = local depth of flow; and b and b, are coefficients whose 
values according to the Manning formula (Eq. 19) are b = 1.0, and 
b^ “ 5/3. Sium (12) investigated the applicability of*^Eq. (19) for 
synthesizing the flow distribution in several natural rivers in- 
cluding the Missouri and Mississippi Rivers. That study indicated 
that the coefficients b and b, in Eq. (19) varied with the ratio 
of width to average depth, B/dT For cross-sections in straight 
river reaches, approximate values of b^ and b^^ are 


bo = 1-0, b^ = 5/3; 50 < B/d < 70 
bo = 0*92, bi = 7/4; B/d > 70 


( 20 ) 


For cross-sections in bends of sinuous and meandering reaches, 
approximate values of b vary between 0.95 and 0.80, and b, 
between 2.48 and 1.78, is B/d increases from 50 to 100. Usually 
the value of b^ must be adjusted to make p = 1 when z = B. 

Fig. 4 shows the variation of measured q/q with d/d for three 
cross-sections in the Missouri Riven- together with the theoretical 
curve corresponding to q/q = (d/d)^' . A comparison of the measured 

and synthesized transverse distribution of unit discharge at one 
of the cross-sections (Mile 410.20) is also shown in Fig. 4. The 
data shown in Fig. 4 correspond to a Missouri River discharge of 
61,300 cfs (1736 cu.m. sec). Synthesized transverse flow distribu- 
tions at Mile 410.20 corresponding to a river flow rate of 35,000 
cfs (991. cu.m/sec) are shown in Fig. 5. To obtain the transverse 
distribution curve for depth d at the discharge of 35,000 cfs (991 
cu.m/sec), the water surface level was dropped by 3.8 ft (1.16 m) 
from that for 61,300 cfs (1736 cu.m/sec). This drop in water 
surface elevation, corresponding to reduction in discharge from 
61,300 cfs (1736 cu.m/sec) to 35,000 cfs (991 cu.m/sec), was 
determined from an estimated stage versus discharge curve for the 
Missouri River at Mile 411. The procedure thus assumes that there 
is no significant change in the cross-sectional shape as the river 
discharge and depth are reduced, and also that there is a unique 
stage-discharge relationship. For alluvial channels, these 
assumptions are not strictly true. Changes in discharge can cause 
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different amounts of aggradation and/or degradation of the bed in 
different parts of the cross-section. Also, seasonal shifts in the 
stage discharge relation may occur that evidently are associated 
with temperature-related effects on the bed configuration, among 
other factors. Unfortunately, alluvial channel phenomena are not 
sufficiently well understood to quatitatively evaluate or signi- 
ficantly improve on these assumptions. 


APPLICATIONS OF THE MODEL 

The solutions of the thermal plume model, Eqs. (9) and (10), were 
used to predict plume patterns in the Missouri River downstream 
from the Fort Calhoun Power Station of the Omaha Public Power 
Station District near Blair, Nebraska, and the Cooper Nuclear 
Station of the Nebraska Public Power District near Brownville, 
Nebraska. Fig. 6 shows the measured (5) and the predicted (7) tem^ 
perature-r ise isotherms for the Fort Calhoun Station corresponding 
to two Missouri River discharges. A comparison of the measured 
(6) and predicted values of AT /6 for the Cooper Nuclear Station 
thermal plumes on four different dates is shown in Fig. 7. Back- 
ground river flow and plant-discharge data corresponding to these 
plumes are given in Table 1. The results shown in Figs. 6 and 7 
indicate that the model predicts thermal plumes in the Missouri 
River that are in quite good agreement with the measured plumes. 

One major capability of the present model is that it permits 
the inclusion of the effect of the cross-channel flow induced by 
channel curvature. Because of this, the model can predict expan- 
sions and contractions of the thermal plume as the river flow 
becomes weighted towards one bank or the other. This feature is 
illustrated in Fig. 8 which shows the predicted thermal plume 
pattern near the Fort Calhoun Power Station for a Missouri River 
discharge of 14,000 cfs (396 cu.m/sec). The computations for this 
case were performed for a 12-raile stretch of the river which was 
divided into ten subreaches, as shown in the insert map in Fig. 8. 


CONCLUSIONS 

A computational model, based on a diffusion equation, has been 
developed for predicting the isotherm pattern for shore-attached 
thermal plumes in rivers. Comparison of the model results with 
field measurements indicates reasonaoly good agreement. Further 
refinements of the model, especially in defining the initial dilu- 
tion factor, need to be made based on laboratory investigations. 
The model can be used to assess the environmental impact of power 
plant thermal discharges in natural rivers. 
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NOTATIONS 


The following symbols are used in this paper: 


A 


a 


3 



area of river flow cross section; 
initial dilution factor; 
width of river flow cross section; 
coefficient used in Eq. 19; 
exponent used in Eq. 19; 
diffusion factor; 
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local depth of flow; 

average depth of river flow cross section; 
overall transverse mixing coefficient; 
acceleration due to gravity; 

Manning coefficient; 

fraction of river flow occupied by plume; 

normalized cumulative discharge; 

volumetric discharge rate of thermal effluent; 

river discharge 

discharge per unit width; 

width-averaged discharge per unit width; 

cumulative discharge from near bank 

radius of curvature of river bend; 

energy gradient; 

depth-averaged local velocity of flow; 

cross-sectional average velocity of flow; 

shear velocity; 

longitudinal coordinate; 

position of virtual source; 

transverse space coordinate; 

excess temperature within plume; 

centerline or bank excess temperature within plume; 
dimensionless transverse mixing coefficient; 
initial excess temperature; 
excess temperature after initial dilution; 
variance of the normal probability density functions, 
Eq. 4 . 
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TABLE I . 

BACKGROUND DATA CORRESPONDING TO THE THERMAL PLUMES SHOWN IN FIGURES 6 AND 7 


Power 

Station 

Date 

River-, 

cfs 

Flow 
/sec ) 

Plant 

cfs 

Discharge 
(m^/sec ) 

Initial Excess 
Temperature 

Op (O^J 

Initial 

Dilution 

Factor 

Ft. Calhoun 

5/3/74 

33,000 

(935) 

801 

(22.7) 

18 (10.0) 

2 

Ft. Calhoun 

12/2/74 

18,500 

(524) 

801 

(22.7) 

16.2 (9.0) 

2 

Cooper 

5/23/74 

48,200 

(1365) 

725 

(20.5) 

7.8 (4.3) 

1 

Cooper 

9/27/74 

37,100 

(1051) 

1352 

(38.3) 

8.6 (4.8) 

1 

Copper 

6/30/75 

57,600 

(1631) 

880 

(24.9) 

10.1 (5.6) 

1 

Cooper 

9/18/75 

66,600 

(1886) 

1088 

(30.8) 

10.5 (5.8) 

1 
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RATIO, d/d 

Figure 4. Transverse Distribution of Dept: 
the Missouri River, and the Rel 
Ship. 
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distance from NEBRASKA BANK, Feet 



15 20 25 30 35 4 0 45 50 


Measured (May 3,i974) 



15 20 25 30 35 


40 45 


Computed 



10 15 20 25 30 35 40 45 50 55 


Measured ( Dece mbe r 2, 1974) 



10 15 20 25 30 35 40 45 50 55 

DISTANCE downstream, Feet x lOO 


ire 6. Comparison of Measured and Predicted Temperature-Rise Isotherms for the 
Missouri River Thermal Plumes Near the Fort Calhoun Power Station, 

Ft. Calhoun. Nebraska. 
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